Photoassociation spectroscopy of cold alkaline earth atoms near the intercombination line 
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The properties of photoassociation (PA) spectra near the intercombination line (the weak transition between 
1 So and 3 Pi states) of group II atoms are theoretically investigated. As an example we have carried out a 
calculation for Calcium atoms colliding at ultra low temperatures of 1 mK, 1 /iK, and 1 nK. Unlike in most 
current photoassociation spectroscopy the Doppler effect can significantly affect the shape of the investigated 
lines. Spectra are obtained using Ca-Ca and Ca-Ca* short-range ab initio potentials and long-range van der 
Waals and resonance dipole potentials. The similar van der Waals coefficients of ground 1 So+ 1 Sq and excited 
1 So+ 3 -Pi states cause the PA to differ greatly from those of strong, allowed transitions with resonant dipole 
interactions. The density of spectral lines is lower, the Condon points are at relatively short range, and the 
reflection approximation for the Franck-Condon factors is not applicable, and the spontaneous decay to bound 
ground-state molecules is efficient. Finally, the possibility of efficient production of cold molecules is discussed. 

PACS numbers: 34.50.Rk, 34.10,+x, 32.80.Pj 



I. INTRODUCTION 

There is a growing interest in the properties of cold 
alkaline-earth atoms. One of the main reasons for this inter- 
est is a possible construction of optical clocks, whose pre- 
cision might exceed that of the current atomic standard of 
time (J H. 

In particular, optical clocks based on an intercom- 
bination transition of alkaline-earth atoms are seen as good 
candidates for the next time standard j3l 001 ■ Increased accu- 
racy of time standards is, for example, desired for the search 
for a time-dependent variation of fundamental constants in 
atomic experiments, which thereby would verify claims based 
on astrophysical data|6|. 

The recent observation of Bose-Einstein condensation 
(BEC) in an Ytterbium gas@] raises hopes for alkaline earth 
atoms, which have similar electronic structure. Such achieve- 
ment would allow a study of cold gasses over a wide range 
of temperatures. Milli- and microkelvin temperatures are 
reached by Doppler cooling on the 1 So— 1 Pi resonance and 
the 1 So— 3 Pi intercombination line, respectively. Nanokelvin 
temperatures are typical for Bose condensates and are reached 
by evaporative cooling. 

Another important reason for interest in alkaline-earth 
atoms is the absence of a nuclear spin in some isotopes. This 
offers an unique opportunity towards a fundamental study of 
Doppler cooling. Moreover, isotopes with zero and nonzero 
nuclear spin allow a comparison of Doppler and sub-Doppler 
cooling 01- The description of atom-atom interactions is 
much simpler for a nuclear spin-less system. In fact, the basic 
theory of photoassociation in strong laser fields fTol ITHl might 
be easier to confirm in alkaline-earth gases than in alkali-metal 
gases, where to date most of the research has been done. 

Scattering of atoms in ground and excited states lead to a 
collisional frequency shift that contributes to the error bud- 
get of an optical clock. A Bose condensate crucially depends 
on atom-atom collisions. Photoassociative (PA) spectroscopy 
fl2l. IT3I fl4l fiH [Trjll is one of the most powerful tools to char- 
acterize these scattering processes. It was developed after the 



success of laser cooling of neutral atoms in the 1980's fl7ll . 

We will focus our investigation on PA spectroscopy of ultra- 
cold alkaline-earth atoms. In the presence of laser light, two 
colliding ground-state atoms, labelled by the scattering state 
"g", absorb a photon forming an excited molecular bound 
state "e" flrj|| . This process is called photoassociation. The 
excited state decays to product states "p" leading to detectable 
loss of atoms from an atomic trap. The variation of the atom 
loss as a function of laser frequency gives the photoassocia- 
tion spectrum. The shape of photoassociation lines not only 
depends on the properties of the colliding atoms but also on 
the temperature and other conditions in a trap (3- 

Phoatoassociation spectra close to the resonance of the 
1 Sq- 1 Pi transition in alkaline-earth atoms were theoreti cally 
analyzed by Machholm et a l iTM Eoll and others El . 
Recently, Degenhardt et al. 12311 measured the photoassocia- 
tion spectra of cold calcium atoms near this transition at mK 
temperatures. Takahashi et al. ll24ll used the photoassociation 
spectroscopy to determined the scattering length of 174 Yb. 

This paper analyzes properties of photoassociation spectra 
near the intercombination line, i.e. laser frequencies close to 
the 1 S'o- 3 Pi transition. A dipole transition between pure sin- 
glet and triplet states is forbidden. However, alkaline-earth 
atom states labelled 3 Pi are not pure triplet states and have 
a small singlet component. This component mostly comes 
from mixing with the nearby 1 Pi state and give rise to a 
weak dipole transition between the 1 Sq and 3 Pi state. Con- 
sequently, the 1 So- 3 Pi atomic line has a small natural width. 
As an example we have carried out calculations for calcium. 

We describe the shape of photoassociation lines with very- 
small natural width and weak laser radiation. The Doppler 
effect as well as the photon recoil must be taken into ac- 
count. This is in sharp contrast with the usual treatment of 
PA IH8II in which these two effects are neglected. Secondly, 
we discuss possible patterns of vibrational levels in pho- 
toassociation spectra near the intercombination line. Close- 
coupled ro-vibrational bound states are obtained using an 
interaction Hamiltonian, which is based on our electronic- 
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FIG. 1: Energy diagram of the photoassociation process. The zero 
energy corresponds to separated atoms in the electronic ground state 
with zero kinetic energy; Ea is the energy of an isolated atom in the 
excited state; E e is the energy of the excited molecular bound state 
|e); Tiuj is the photon energy; e r is the kinetic energy of the relative 
motion of colliding atoms; A e is the binding energy of the bound 
state |e); A is the detuning of the photon energy from the isolated 
atom excitation energy Ea- 



structure potentials ll25ll and recently calculated dispersion 
coefficients The Hamiltonian also includes coupling be- 
tween 1 Pi and metastable 3 Po,i.2 states. The interatomic po- 
tential between two ground state calcium atoms is relatively 
well known[27]. It is shown that when the interaction in the 
ground and excited state is similar the reflection approxima- 
tion Jiol fl4l Eijll cannot be applied to calculate the intensi- 
ties of photoassociations lines. This is unlike photoassocia- 
tion spectra near strongly allowed transitions, where the in- 
teraction in the ground and excited states differs significantly 
and the reflection approximation is well satisfied. Finally, we 
show that efficient production of cold molecules in the ground 
electronic state using photoassociation should be feasible. 



II. SHAPE OF THE PHOTOASSOCIATION LINE 

The photoassociation process occurs in a thermal cloud of 
cold atoms at temperature T interacting with weak laser radi- 
ation. After absorption of a photon of frequency cu, two atoms 
form an excited molecular bound state |e) with energy E e . 
The photoassociation process according to standard descrip- 
tions is most efficient when the photon energy and the kinetic 
energy of the relative motion of colliding atoms, e r , match the 
energy of the excited bound state, that is hco + e r = E e . 

A schematic of energies in the PA process is shown in 
Fig-E The binding energy A e of bound state |e) is given by 
A e = E e — Ea, where Ea is the energy of an isolated atom 
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FIG. 2: Schematic of the collision before the photoassociation pro- 
cess in the laboratory frame and in the center of mass coordinates. 
pi and p2 are momenta of the colliding atoms; p r is the relative mo- 
mentum of the colliding atoms; p c is the center-of-mass momentum 
of the two atoms in the laboratory frame; Tiuj is the photon energy 
in the laboratory coordinates; Tiuj + ed(p~c) is the Doppler shifted 
photon energy in the moving frame of center of mass coordinates. 



in the excited state and the zero energy corresponds to sep- 
arated atoms in the electronic ground state with zero kinetic 
energy. Furthermore, the detuning A of the photon is defined 
by A = Tiuj — Ea, where Tiuj is the photon energy. The reso- 
nance condition for the PA process in terms of detunings reads 
A + e r - A e = 0. 

The description of photoassociation for very narrow lines 
requires us to include two new effects. These effects are 
the Doppler shift and the photon recoil. Figure [2] shows a 
schematic of the collision before the PA process. In the lab- 
oratory frame the two colliding atoms, each with mass m, 
have momentum pi and p-2, respectively. In the coordinate 
frame, which moves along with the center of mass, the rela- 
tive momentum of the colliding atoms is p r and kinetic en- 
ergy of relative motion e r (p r ) = p^./2/j,, where [i = m/2 
is the reduced mass of the colliding atoms. The photon en- 
ergy hu> + £d(Pc) in the moving frame is shifted with respect 
to its energy Tiuj in the laboratory frame. To good approx- 
imation the Doppler shift £d(Pc) = — Tifcias ■ Pc/M is pro- 
portional to p c , where p c is the center-of-mass momentum of 
the two atoms in the laboratory frame. The total mass of the 
system is M = 2m, fci as is the wave vector of the laser ra- 
diation with absolute value fci as = ui/c, and c is the speed of 
light. After photoassociation the excited molecule gains the 
momentum Tik\ as of the absorbed photon and, therefore, has 
a translational kinetic energy of £/ re c,moi = ^ 2 kf as /(2M) in 
the moving frame defined before the absorbtion of the pho- 
ton. Consequently, photoassociation is most efficient when 
Tiuj + £d(Pc) + £ r (pr) = E e + E ICC , mo \ or in terms of detun- 



ings A + £ D {fc) + £r{fr) ~ A e - E T , 



,mol 
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The excited molecular states created in the PA process ei- 
ther decay back to the ground state or can be further excited 
to ionizing states. In the former case the product states escape 
from the trap and give rise to trap loss. Ions are detected in the 
latter case. In this paper we will only model trap loss. The loss 
mechanisms are characterized by a rate coefficient K(A,T), 



3 



which describes the efficiency of the process for a given laser 
detuning, intensity /, and atom temperature. For clarity we 
omit / as argument in the rate coefficient. The rate coeffi- 
cient K(A, T) is linear in the weak laser field intensity /. For 
hi ghe r intensities the center of the PA line shifts linearly with 
I I10tf29tl30ll and the width of the line increases. The descrip- 
tion of the PA line shape presented in this work, for simplicity, 
omits this shift as it can be neglected in the weak laser field 
regime. 

The photoassociation trap-loss rate coefficient involves a 
thermal average of the rate coefficient for a given pair of mo- 
menta pi and p2 in the laboratory frame. The momentum 
distribution of both atoms is Maxwellian with temperature T. 
In practice we average over the momenta p c and p r , denoted 
by (...). These momenta have a Maxwell-Boltzmann distri- 
bution / M (p c ) = {^/np M )~ 3 exp(-p 2 c /p 2 M ) and f»(p r ) = 
(y/irp^)' 3 exp(— Pr/pft) with the most probable momentum 
Pm = ^/WzbTM and p^ = ^/2fcgT/I, respectively, where 
ks is the Boltzmann constant. Then, we have 

K(A,T) = (IC(A,p c ,f r )) = 

J d 3 p c f M (Pc) J d 3 p r f^Pr) K.(A,p c ,p r ) , (1) 

where JC(A,p c ,p r ) describes the trap loss from a collision 
with momenta p c and p r . 

The trap-loss coefficient for p c and p r is equal to 

K(A,p c ,p r ) = ^2 v rT^\S pg {A,p c ,p' r ; e)| 2 , (2) 

where k r is the relative wavenumber defined by e r = 
h 2 k 2 /(2/j,) and v r = hk r /fi is the relative speed of the col- 
liding atoms. The quantity \S pg (A,p c ,p" r ; e)| 2 is the transi- 
tion probability from an initial ground state, g, to all product 
states, p, 16611 through an intermediate excited bound state e. 
The indices g and e, summed over in Eq. l|2}, represent quan- 
tum numbers that describe the initial and intermediate state, 
respectively. Each initial state g is labeled by the total angular 
momentum quantum number J g , its projection quantum num- 
ber M g , and the total parity p g , while the intermediate rovi- 
brational levels e are labeled by vibrational quantum number 
v, total angular momentum J e , projection M e and parity p e . 

The transition probability from an initial state to the product 
states can be described by a generalized resonance formula 

\S pg (A,tf c ,p r ;e)\ 2 = 

(A + e D {f c ) + e r {pr) - A e - ^rccanol) 2 + (r e /2) 2 ' 

(3) 

where the total width of the excited bound state e, T e = 
r e ,„ at + r eidis + J2 g T eg(£r), is the sum of its natural ra- 
diative width r e nat , the contribution T e di S from predissoci- 
ation, and the stimulated widths Y eg {e r ) caused by the laser 
coupling between the excited and ground states. The width 
T pe describes the decay into product states and is assumed to 



be equal to the sum of the natural radiative and predissociation 
width of the bound state. 

The stimulated width is proportional to the light intensity 
and is calculated from Fermi's golden rule llOlBlll 

T ea {e r ) = 27r|(* e («, J e M e p e )\V\ &s \*+{e r ,J g MgPg))\ 2 , 

(4) 

where (^(w, J e M e p e )) is the unit-normalized excited bound 
state and \^^(e r) J g M g p g )) is the energy normalized scat- 
tering ground state. The operator V\ as describes the coupling 
between the ground and excited state by laser light. Details of 
the close-coupled equations l32ll that are solved to calculate 
the bound and scattering states are described in appendix lAl 
To describe the atom-photon interaction during a collision we 
adopt the treatment developed by Napolitano et al. l33ll . The 
matrix elements of V[ as are given in appendixlBl 

To highlight properties of the photoassociation lines we re- 
duce the thermal average in Eq. Q to the 2D integral 

(^j^r\Spg{A,p c ,p r ;e)\ 2 ^ = 

+ OO OO 

/ [dxxe-* 2 C(A,y,x 2 ), (5) 

hQr J J 

-00 

where y = -k las ■ p c /(k ias p M ) and x = p r /p P are di- 
mensionless variables similar to those used in the description 
of pressure and Doppler broadened spectral lines l34l l35ll . 

Qt = (2n^k B T/h 2 ) 3 / 2 , and 

£(A,y,x 2 ) = 

r pe r eg (x 2 A T ) 

(A + yA D + x 2 A T - Ae - £ rec ,moi) 2 + (r e /2) 2 - 

(6) 

The quantities At = k B T and Ad = hk\ as ^2kBT/M are 
the thermal and Doppler width, respectively. 

Three limiting cases of the lineshape, Eq. are of in- 
terest. The shape of the line is Lorentzian when T e is much 
bigger than Ad, At, and T eg (e r ). The denominator of C can 
then be pulled out of the integrals leading to a Lorentzian pro- 
file with a full-width at half maximum (FWHM) equal to T e . 
Such a lineshape can be expected for strongly allowed tran- 
sitions at ultra-low temperatures on the order of nanokelvins, 
such as exist in Bose condensates I36ll37ll l67ll . 

In a second limiting case the shape of the line is a "cut- 
off exponential". This profile can be obtained when At is 
much larger than Ac and T e and the energy dependence 
of T eg (e r ) is neglected (compare discussion in Ref. fiHl 

where the energy dependence T eg (e r ) ~ s l r +1 ^ 2 describ- 
ing the Wigner threshold law behavior is considered). The 
Lorentzian C can be replaced by a delta-Dirac function with 
argument A — A e — E Iec>mo \ + x 2 At and the integrals can 
be solved analytically. In fact, the profile is proportional to 
6{— A) cxp(A/At), where 6{z) is the Heaviside step func- 
tion: 6{z) = 1 for z > 1 and 9{z) = for z < 1. The full 
width at 1/e of the exponential lineshape equals the thermal 
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width At- This lineshape is most easily observed at magneto- 
optical trapping temperatures on the order of a millikelvin. 

Finally, for the unusual situation of extremely weak tran- 
sitions at nanokelvin temperatures one could try to achieve 
conditions in which is much bigger than Ay and T e . In 
such a case the Lorentzian can again be replaced by a delta 
function but now with an argument that only depends on y 
and A. The resulting lineshape is a Gaussian with half width 
at 1 Je of the maximum equal to Ad . 

In the usual treatment of the PA line shape Doppler broad- 
ening is neglected. To find conditions for which this ap- 
proximation breaks down, the relative importance of Doppler 
and thermal effects must be determined. It is easy to show 
that Ad = At at temperature T = Tn, where Tr = 
h 2 kf as /(mkB) is the atomic recoil temperature. At temper- 
atures T > Tr thermal broadening dominates, while for 
T < Tn Doppler broadening can determine the shape of the 
line. In fact, this requirement is not sufficient. It is also neces- 
sary to assume that the Doppler width is comparable or bigger 
than T e . 

The influence of the photon recoil on the PA spectra is 
to very good approximation described by an uniform shift, 
£rec,moi, of all PA lines. It should be noted that the pho- 
ton recoil energy of the two-atom molecule is two times 
smaller than the photon recoil energy of an isolated atom 
E rcc = h 2 kf as /(2m). This indicates that for molecular bound 
states close to atomic levels there should exist a transition re- 
gion from recoiling as a molecule to recoiling as an atom. Our 
theory does not treat this effect and should only be applied to 
molecular states with a binding energy that is much bigger 
than the photon recoil energy. 



III. INTERATOMIC HAMILTONIAN 
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FIG. 3: Adiabatic potentials of two calcium atoms as a function of 
interatomic separation. The inset shows a blowup of the region near 
the 3 P + 1 S limits. At short range the potentials are labelled by 
their Hund's case (a) symmetry and in the inset by a Hund's case (c) 
symmetry label. Ungerade excited states, accessible by an optically 
allowed transition from the ground X 1 E ! j" potential, are shown. 



Properties of photoassociation spectra are governed by the 
interactions between colliding atoms. The interaction Hamil- 
tonian is described in appendix 1X1 and is similar to that dis- 
cussed by Mies et al. 13911 for the electronic structure and spec- 
troscopy of Hg2. This Hamiltonian includes non-relativistic 
Born-Oppenheimer potentials, the spin-orbit splitting of the 
3 -Po,i,2 atomic states, relativistic coupling between 1 P and 
3 P states, and a term that incorporates the rotation of the two 
atom system. 

The best possible Born-Oppenheimer potentials are used. 
Ab initio calculations, with the exception of a few simple 
cases, do not give sufficiently accurate predictions of abso- 
lute positions of molecular bound states. Therefore theoreti- 
cal potentials are used as initial guesses, and modified at short 
interatomic separation to reproduce experimental binding en- 
ergies. Unfortunately there is no experimental data on binding 
energies near the 3 Pi + 1 Sq dissociation limit for the calcium 
molecule. As the rovibrational structure is not known experi- 
mentally, we can only map out possible spectra by varying the 
short-range part of the potential. 

The Born-Oppenheimer potentials have Hund's case (a) 
symmetry, 2S+1 |A| (T , where S is the total electron spin, A is 
the projection of the total electron orbital angular momentum 



along the interatomic axis, and a = g/u describes the gerade 
or ungerade symmetry of a state. Near the 3 Pi + 1 Sq limit, 
where the atomic spin-orbit interaction is much bigger than 
the Born-Oppenheimerpotentials, adiabatic potentials are bet- 
ter described by Hund's case (c) symmetry, |0| J, where f2 is 
the projection of the total electron angular momentum along 
the interatomic axis and for = 0, the label ± describes the 
symmetry under a reflection of the electronic wavefunction. 
The adiabatic Hund's case (c) potentials are obtained by si- 
multaneous diagonalization of the Born-Oppenheimer poten- 
tials and the spin-orbit coupling. 

Relevant data about the interaction potentials between 
alkaline-earth atoms in the ground and excited states have 
been compiled by Kotochigova and Julienne 12^1 . The Ca2 
potentials l25ll are shown in Fig. [3] Czuchaj et al. l40ll have 
published a similar set of potentials. It is convenient to use 
atomic units. The atomic unit of length 1 au = 1 «o = 
0.0529177 nm is equal to a Bohr radius ao = Ti 2 / (m e e 2 ); the 
atomic unit of energy 1 au = E h = 4.3597482 x 10~ 18 J is 
equal to the Hartree energy Eh = h 2 / (m e ag). Here m e and 
e are the electron mass and charge, respectively. 

The A 1 ^^ electronic ground state potential of Ca2 has 
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been determined by Allard et al. fiFh . Parameters of the po- 
tential are listed in Table 1 of this reference. For large inter- 
atomic separations the potential asymptotically approaches a 
van der Waals potential with C 6 (X 1 E+) = 2081.18 au. The 
atomic unit of C 6 is 1 au = E h a% = 0.957342 x 10~ 79 Jm 6 . 
The scattering length for this potential is a scat = 389.8 au. 

The exited-state interaction potentials between calcium 
atoms dissociating to the 1 Pi + 1 5o and 3 Po,i.2 + 1 «5o limits for 
short-range interatomic separations have been modelled using 
the adiabatic 1,3 E+ U and 1,3 n g u potentials. These adiabatic 
potentials are determined on the basis of ab initio calculations 
of Ref . l25ll and smoothly connected to their asymptotic func- 
tional form. The form is Ce / r 6 for triplet states and C3 /r 3 
for singlet states, respectively. Triplet potentials dissociate to 
the 3 P + 1 S limits, while singlet potentials dissociate to the 
X P + 1 S limit. 

In the model describing interaction of two calcium atoms 
we do not include potentials correlating to 1 Z?2 + 1 So lEUl and 
3 -01,2,3 + 1 Sq dissociation limits. These potentials could give 
rise to molecular bound states near the 3 Pi + 1 So limit, but the 
sparse density of states makes it unlikely that such level occur 
in a small energy interval close to the 3 Pi + x 5q dissociation 
limit. 

Dispersion coefficients for two ground state atoms are 
well known for many atomic species 1271 l42l 14311 . For ex- 
cited atoms, however, there is little data. For a calcium 
atom in the ground 1 S state interacting with another calcium 
atom in the excited 3 P state the long-range dispersion co- 
efficients have been recently calculated by Derevianko and 
Porsev l26Tl . The attractive Hund's case (c) |0|^, u poten- 
tials correlating to the 3 Pi + 1 So limit have Ce(0^,) = 

2462 au and C 6 (l s/tt ) = 2593 au 1 26]. For this Pa- 
per, however, dispersion coefficients are needed for Hund's 
case (a) 2S+1 |A| 9 / U Born-Oppenheimer potentials. Follow- 
ing Ref. |H we have C 6 (0+J = C 6 ( 3 IT s/u ), C 6 (l g/u ) = 
[C 6 ( 3 E S/U ) + C 6 ( 3 n 9/u )]/2, and, therefore, C 6 ( 3 E 9/tl ) = 
2724 au, C 6 ( 3 II 9/u ) = 2462 au. Derevianko and PorsevlH 
found that the splitting between the long-range E and II po- 
tentials is small. This small anisotropy for the quasi-two- 
electron atom is qualitatively different from analytical predic- 
tions for single-electron excited alkali-metal atoms interacting 
with rare-gas atoms El Efj ll . Older data l46ll does not agree 
with the results of Ref. 12611 . 

The lifetime of the atomic 1 Pi state determines the dipole- 
dipole interaction coefficients, C3, of the singlet poten- 
tials. In fact, C3OIAW = COlAWCf 5 where C*f } = 
(3/4)r A ( 1 P 1 )[A A ( 1 Pi)/(2 7 r)] 3 , T^( 1 P 1 ) = n/T^Px) is 
the natural width of the excited 1 P\ state, and Xa is the 
wavelength of the corresponding radiation. The coefficient 
COIAW is defined by C(%) = +2, CO^u) = -2, 
COn,,) = -L and C^ILj) = +1. For calcium r A ( 1 Pi) = 
4.59 nsJHil. 

Retardation effects EH l49ll do not change our main conclu- 
sion. The calcium bound states of interest are mostly con- 
fined to interatomic separations that are small compared to 
the wavelength Aa( 1 Pi)- Under such circumstances retarda- 
tion effects can be neglected. Moreover, this implies that pho- 
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FIG. 4: An example bound-state structure and corresponding pho- 
toassociation spectrum of calcium near the intercombination line. 
The detuning A is defined relative to the 3 Pi + 1 So limit. The panel 
on the left shows the 0„ (dark line) and l u (light line) Hund's case 
(c) adiabatic potentials as a function of interatomic separation r. The 
bound states of these potentials are indicated by the thick and thin 
solid horizontal lines, respectively. The panel on the right shows the 
corresponding photoassociation spectrum under typical conditions. 
Each peak corresponds to a bound state in the left panel as indicated 
by the dashed lines. Finally, the short dotted lines on the left side of 
the graph indicate "bin edges" of the 0^ potential. 

toassociation of two ground state atoms can only excite the 
ungerade states. Therefore, Fig. [5] only shows ungerade ex- 
cited potential curves. 

IV. RESULTS 

Photoassociation spectra near the 3 Pi + 1 5'o limit are ex- 
pected to be weak because the atomic transition dipole to the 
intercombination line is nearly forbidden. As we will show, 
such spectra can most easily be measured at ultracold temper- 
atures on the order of /iK and below. At these temperatures 
only s-wave collisions will contribute to the spectrum. Conse- 
quently, contributing transitions are between the ground ger- 
ade scattering state, |^+(£ r , J g M g p g )), with total molecular 
angular momentum J g = l g = and parity p g = 1, where l g 
is the rotational angular momentum between the atoms, and 
excited ungerade bound states | v I / e(i', JeM e p e )) with J e = 1 
and p e = — 1. For total molecular angular momentum J g = 
and J e = 1 there are one and five coupled channels, respec- 
tively (See AppendixlAl). The excited-state channels for odd 
J e and negative total parity p e do not include the 3 Po atomic 
state and, therefore, the predissociation width r e di s in Eq. [4] 
is zero l6Hl . 

Figure |4] shows an example of a bound state structure and 
the corresponding spectrum. The PA spectrum is a reflection 
of the rovibrational structure of excited molecules. In this case 
lines are assigned to the Hund's case (c) 0+ or l u symmetry. 

Details of the Born-Oppenheimer potential at short inter- 
atomic separation are insufficiently known for a quantitative 
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prediction of the bound-state locations. Therefore, we have 
modified the short range of the 3 n u potential to demonstrate 
how a spectrum can change. The 0+ and l u potentials corre- 
late to the 3 n u and 3 £+ potentials at short range, respectively. 
Consequently, by changing the 3 II [I potential we can change 
the location of the 0+ bound states while leaving the position 
of l u bound states virtually unchanged. It is convenient to de- 
fine "bins" of the 0+ potential as energy intervals with edges 
marked by the 0+ energy levels calculated with a 3 II U poten- 
tial, such that the last bound state is exactly on the 3 Pi + 1 Sq 
threshold. As the short-range 3 H U potential is changed, there 
is always exactly one 0+ bound state in each bin. The "bin 
edges" for 0+ states are shown in Fig.|4] 

Figure |5] shows PA spectra for two 3 n„ potentials in or- 
der to illustrate the limiting cases of overlapping and non- 
overlapping 0+ and l u bands. The PA spectra are calculated 
using the lineshape formula derived in the previous section 
and take into account Doppler broadening. In both spectra the 
last five vibrational levels of 0+ and l u symmetry are shown. 
The PA rate coefficient increases by nine orders of magnitude 
near each vibrational level. Observable trap loss is on the or- 
der of the 10~ 12 cm 3 s _1 . In principle, ion detection allows 
measurement of weaker PA lines than is possible with trap 
loss. On the frequency scale of the figure the location of the 
l u vibrational lines for the two potentials is almost the same. 

A spectrum in which lines with 0+ and l u symmetry are 
far apart from one another is shown in Fig. [2a). In this case 
the projection O e of the electronic angular momentum j e on 
the intermolecular axis is a good quantum number and bound 
states can be labelled by the Hund's case (c) coupling scheme. 
This kind of spectrum has been observed in preliminary exper- 
iments on Strontium near the 1 So— 3 Pi line by Ido and Katori 

The rather small anisotropy of the long-range dispersion in- 
teraction of the 0+ and l u potentials (Ce(0+) w Cq(1 u ) to 
within 5%) can lead to a near coincidence of 0+ and l u lev- 
els over a range of v. Such a case is shown in Fig. 0b). For 
closely spaced doublets the projection fl e is not a good quan- 
tum number, and bound states should rather be labelled by the 
rotational angular momentum l e as in the Hund's case (e) cou- 
pling scheme. It can be clearly seen that in each doublet there 
is a strong and weak line corresponding to l e = and l e = 2, 
respectively. 

We have also studied the change in coupling scheme from 
Hund's case (c) to e for bound states very close to the molec- 
ular thresholds. By changing the shape of the 3 n„ poten- 
tial, the last bound state, initially attributed to 0+ symme- 
try, smoothly approaches threshold. Simultaneously, the wave 
function smoothly changes its character as well. For a binding 
energy A e /h = —0.026 GHz 90% of the wave function has 
Q e = character, while for A e /h = -0.003 GHz 90% of 
the wave function can be attributed to l e = 0. In other words, 
Hund's case (c) holds for binding energies larger than 0.026 
GHz. 

In Fig. 13 and in the remainder of this paper the natural 
linewidth r e nat is approximated by 2Ta( 3 Pi) — 0.663 kHz. 
The natural linewidths calculated from the theory described 
in the appendices shows that the width is less than about four 
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FIG. 5: Photoassociation spectra of calcium near the intercombina- 
tion line for two 3 I1„ potentials. Panel (a) is for the potential used in 
Fig-E] Panel (b) is for a potential where the and l u bands nearly 
overlap. The laser intensity is 1 W/cm 2 and the temperature of the 
gas is T — 1 /iK. The and l u vibrational assignment is shown in 
Panel A. The last vibrational level is labeled by "1". 



times r^( 3 Pi) for the detunings shown in Fig. [3] We believe 
that our simple model can not quantitatively describe the natu- 
ral widths. The treatment of the coupling between 1 P and 3 P 
states is not sufficiently accurate. For simplicity we assume 
the same natural linewidth for all lines. 

Figure displays the collision-energy dependence of the 
stimulated width r eg (£ r ) for the 0+ lines shown in Fig. 0a). 
The ground-state potential has a scattering length of 389.8 ao 
and the laser intensity is 1 W/cm 2 . The width rapidly in- 
creases for bound states closer to the threshold as the over- 
lap of the bound and scattering wavefunction grows. More- 
over, for collision energies smaller than e r /ks = 100 yuK the 
width of all lines is proportional to t/s^, satisfying the Wigner 
threshold law. For a 1 W/cm 2 laser intensity and collision en- 
ergies less than 1 fj,K the stimulated width for four of the vi- 
brational levels is smaller than the natural width. Therefore, 
the lines are unsaturated for most of lines shown in Fig. [5] The 
exception are those lines closest to resonance. 

In order to show the temperature dependence of the shape 
of a line we have chosen line 2 of the 0+ band in Fig. 0a). 
Figure Q shows the lineshape for a temperature of 1 mK, 1 
/iK, and 1 nK with and without Doppler broadening. In the 



7 




10"" 10 

e/k B [K] 

FIG. 6: The stimulated width as a function of the kinetic energy of a 
S + 1 S collision calculated for the five lines of the 0„ band shown 
in panel (a) of Fig. [5] The laser intensity is 1 W/cm 2 . The natural 
width r e , na t/ft = 0.663 kHz is marked by the dashed line. 
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FIG. 7: The lineshape of the 0„ line labelled 2 shown in panel (a) 
of Fig. [5] as a function of laser detuning. Panels (a), (b) and (c) 
show lineshapes at a temperature of 1 mK, 1/jK and 1 nK, respec- 
tively. Profiles that include and do not include Doppler broadening 
are marked by solid and dashed lines, respectively. The laser inten- 
sity is 1 W/cm 2 . The quantity A e is the binding energy of the 0^" 
rovibrational level relative to the So + 1 Pi dissociation limit. 



absence of Doppler broadening Ac = in Eq. |6] For these 
three temperatures the thermal width At/Ii = 20837 kHz, 
20.8 kHz, and 0.021 kHz and the Doppler width A D /h = 
694 kHz, 21.9 kHz, and 0.694 kHz, respectively. 

At a temperature of 1 mK the line in Fig. 0a) has the typical 
"cut-off exponential" shape determined by thermal broaden- 
ing and is only slightly affected by Doppler broadening. The 
width of line is on the order of ten MHz and the peak rate co- 
efficient is 10 -15 cm 3 s _1 . Such a low rate coefficient makes 
trap loss hard to detect in typical ultra-cold-atom experiments 
in a magneto-optical trap. In this case ion detection might be 
a sensitive alternative. 

A 1 /iK atomic-gas temperature is close to the recoil tem- 
perature Tr = 1.11 /iK. Under such conditions both thermal 
and Doppler broadening in the PA lineshape are comparable. 
Figure Qb) demonstrates a significant difference between a 
Doppler broadened profile and one without Doppler broaden- 
ing. The width of the line is on the order of hundred kHz and 
the peak rate coefficient is 10 -12 cm 3 s" 1 . A trap-loss signal 
should be observable for such rate coefficients. 

In Fig. 0c) a PA lineshape for a thermal gas at 1 nK is 
shown. Typically for such low temperatures and sufficiently 
high densities an atomic gas could be Bose condensed(BEC) 
and a Boltzman distribution of atomic momentum should not 
be used. Here, we assume a low enough density that conden- 
sation has not occurred. The lineshape is an ordinary Voigt 
profile, which is determined by Doppler and natural broad- 
ening. The natural width r e nat is 0.663 kHz. The line is 
Lorentzian if Doppler broadening is neglected. The width 
of the line is on the order of kHz and the peak rate coeffi- 
cient is 10~ 10 cm 3 s _1 . Moreover, the molecular recoil en- 
ergy E rcC:mo i/h = 5.775 kHz is significantly bigger than the 
width of the line. 

The width of the lines in Fig. Q varies by four orders of 
magnitude. The peak rate coefficient changes by five orders of 
magnitude. Clearly, for temperatures on the order of 1 yuK and 
below, photoassociation spectra should be observable. More- 
over, Doppler broadening is an important factor and affects 
the shape of the lines significantly. 

We have also analyzed the stimulated width of the J lines 
as a function of the scattering length in the 1 Sq + 1 So ground 
state and the binding energy A e of the excited 0+ bound 
states. The scattering length is varied within the range allowed 
by experimental. The scattering length is varied by slight 
modifications of the short-range part of the ground-state po- 
tential. The binding energy of the 0+ bound states are changed 
by modifying the 3 n u potential as discussed above. 

Figure [3] shows the stimulated width at a collision energy 
Er/ks = 1 as a function of a scat and A e . The fig- 
ure shows multiple nearly vertical dark structures, where the 
width is nearly zero. There are two kinds of these structures: 
ones accompanied by a parallel bright feature, where T eg is 
large, and those without. The latter structures are not quite 
vertical for smaller a sca t. The "first kind" of structure occurs 
when a 0+ bound state coincides with a l u bound state. This 
mixing is independent of any ground state scattering property 
and therefore the structures are vertical. The "second kind" of 
dark structures occur when the overlap of the scattering wave- 
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FIG. 8: The stimulated line width of the 0„ band as a function of 
their line position (binding energy) A e and scattering length a scat of 
the ground state potential. The laser intensity is 1 W/cm 2 and the 
kinetic energy of the collision is e r /fcs = 1 AtK. 



function and the excited bound state vanishes. The shape of 
the scattering wavefunction near the outer turning point of the 
excited bound state does not change much when the scattering 
length is on the order of a few hundred an. Therefore, these 
structures are nearly vertical. Near a scat = 300 o,q a departure 
from vertical can be observed. 

For photoassociation near strongly-allowed atomic tran- 
sitions, such as occur in alkali-metal experiments fTM fT3l 
and near the 1 P\ line of the alkaline earths 1 1911 . the overlap 
vanishes for excited bound states with outer turning points 
near the nodes of the ground-state wave function. Their 
detuning A e can be found with the help of the reflection 
approximation (3 HHh which says that the Franck-Condon 
factor is proportional to the square of the ground state wave- 
function at the position rc for which the difference in the ex- 
cited and ground state potentials equals the photon frequency. 

In our case the reflection approximation cannot be ap- 
plied because the asymptotic potentials in the ground and ex- 
cited state are similar. The ground-state wave function has 
nodes at 26.5 au, 32.4 au, and 46.1 au for a scattering length 
o-scat = 389.8 au and a collision energy e r /kB = 1 A*K. 
For the reflection approximation the relation between detun- 
ing and outer turning point of the excited bound states is de- 
termined from a potential that is the sum of the Hund's case 
(c) potential V(0+) and the rotational correction [2 + J e (J e + 
l)]h 2 /(2/w 2 ). The potential V(0+) approaches the 3 n u po- 
tential at large internuclear separation and J e = 1. These 
ground state nodes correspond to outer turning points detun- 
ings of -51.3 GHz, -15.0 GHz, and -1.6 GHz, respectively. 
These detunings do not correspond with dark lines in Fig. [S] 

Appendix Icl gives another perspective of Fig. [8] The data 
are described in terms of a near threshold vibrational quantum 
number instead of the binding energy. Integer values of this 
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FIG. 9: The fraction f ge of molecules in rovibrational state v g ,J g = 
of the EJ" potential produced by natural decay of rovibrational 
state v e ,J e — 1 of 0„ symmetry. 



quantum number are related to the bins defined in Fig. [4] This 
discussion is not crucial for the main thrust of the paper and, 
therefore, has been placed in an appendix. It, however, gives 
a deeper understanding of the physics involved and is worth 
presenting. 

Finally, we have investigated the possibility of creating cold 
molecules in the ground electronic state via the photoassoci- 
ation process. The similar Cs/r 6 dependence of the ground 
and excited potentials can make such a process more efficient 
than making ground state molecules through excited states in 
which the asymptotic form of the potential is Cs/r 3 . We have 
calculated rates for bound-bound transitions between 0+ ex- 
cited states and 1 E+ ground states. From these rates we have 
calculated the fraction of ground-state molecules, which af- 
ter photoassociation are formed by natural decay. Figure [9] 
presents the fraction f ge of molecules in rovibrational state 
V g , J g = of the potential produced by natural decay of 
rovibrational state v e ,J e = 1 of 0+ symmetry. The fraction 
of a v g , J g = molecule is no more than 15%. For v e > 10 
high Vg vibrational levels of the ground state potential are pro- 
duced, while for v e < 10 low v g levels are produced. Never- 
theless, even for high v e a few percent of the ground-state 
molecules have low v g . The fraction of molecules decaying to 
Jg = 2 vibrational states of the 1 S g potential have also been 
calculated, but are not shown. It was found that the fraction 
of excited molecules decaying to the J g = 2 vibrational states 
is bigger than those decaying to the J g = vibrational states. 
Rotational states with J g = 1 do not exist because of Bose 
symmetrization (See appendixlATi. 

Adding the J g = and 2 fractions f ge of molecules shows 
that most of the spontaneous emission of a v e , J e 0+ rovibra- 
tional level goes to ground molecular bound states. In fact, for 
v e < 19 more than 80% of the decay is to molecular states. 
The only exception is the last vibrational level, v e = 19, cor- 
responding to line 1 in Fig.^Ja), for which 30% of the decay is 
to molecular states. We have also found that the spontaneous 
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decay rate of deeply-bound excited states can be one order of 
magnitude larger than 2Ta{ 3 Pi)/^i- 



V. CONCLUSIONS 

We have developed a description for the shape of photoas- 
sociation lines for weak transitions in which the natural width 
is smaller than the Doppler width. The lineshape theory in- 
cludes Doppler broadening and a photon recoil shift. It was 
shown that the Doppler effect significantly affects the PA line- 
shape when the gas temperature is on the order of the recoil 
temperature and below. 

A model calculation for calcium has been carried out. It 
is an example of photoassociation near the intercombination 
line of alkaline-earth atoms. We find that photoassociation 
spectroscopy should be possible at calcium gas temperatures 
on the order of fiK and below. In addition, it was shown that 
when the long-range potential of the ground and excited state 
are similar, the reflection approximation fiol fl4l l28ll incor- 
rectly estimates the stimulated width and strength of PA lines. 
Finally, we have indicated that excited molecules are very 
likely to decay to vibrational levels of the electronic ground 
state. This can be used as an effective way to produce cold 
molecules. 
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APPENDIX A: CLOSE-COUPLING CALCULATIONS 

The quantum theory of slow-atom collisions I5TII52I1 allows 
a quantitative description of the scattering process and bound 
states. We apply this theory to describe the slow collisions 
and bound states of alkaline-earth atoms. 

Scattering or bound states of two identical alkaline-earth 
atoms without nuclear spin can be described in the basis 

\SLjlJMj;p)= (jlmjmilJMj^SLjmi^Yi^e,^) 

rrij ,mi 

(Al) 

where the \SLjrrij; a) describe the electronic state of the 
molecule and Yi m (9,4>) is a spherical harmonic describing 
the nuclear rotation. The quantity L is the total electron or- 
bital angular momentum, S is the total electron spin angular 
momentum, j = L + S is the total electron angular momen- 
tum, / is the rotational angular momentum, and J = I + j is 
the total angular momentum. The projections of j, I, and J on 
a space-fixed z-axis are rrij, mi, and Mj, respectively. The 



quantity p is the total parity. Gerade (g, a = +1) and unger- 
ade (u, a = — 1) electronic states correspond to total parity 
p = +1 and p = — 1 states, respectively. This is a conse- 
quence of the more general rule for atoms with nonzero nu- 
clear spin p = a(— 1) 2/ , where I is the atomic nuclear spin. 
The total parity restricts the allowed I by p = PaPb ( — 1) > 
where pa and pb are the atomic parities. The atomic parity 
is +1 for the ground 1 S state and —1 for excited 1 P and 3 P 
states. 

The molecular Hamiltonian H = T + Ha + V- ln t + 
Vrot is calculated in the \SLjlJMj;p) basis. Here, T = 
— (h 2 /(2fj,))d 2 /dr 2 is the kinetic energy operator, Ha is 
the atomic Hamiltonian, Vint are the nonrelativistic Born- 
Oppenheimer potentials, and V YO t = fir j (2fj,r 2 ) describes 
the rotational energy. The matrix elements for the kinetic 
and rotational energy are diagonal in this basis. In fact, 
(SLjlJMj;p\T\SLjlJMj;p) = -(h 2 /{2 l j))d 2 /dr 2 and 
{SLjlJMj;p\V Iot \SLjlJMj;p) = K 2 l{l + l)/(2 M r 2 ). 

The matrix elements for the Born-Oppenheimer potentials 
are calculated in two steps. The first step involves transform- 
ing the molecular electronic state \SLjrrij; a) into a body- 
fixed coordinate system. That is into a superposition of 
\SLjQ,; a) states, where fl is the projection of j on the in- 
ternuclear axis. After some algebra the matrix elements are 
given by 

(S'L'j'l'J'M'j;p'\Vi nt \SLjlJMj,p) = 

/ (2Z' + 1)(2Z + 1) 
^, p 5j,,j5 MIjMj y (2J+1)2 

(j'l'VO I fO) {jlSlQ I JO) (S'L'j'n; a\V int \ SLjSl; a) , 

n 

(A2) 

i.e. the operator Vj n t is diagonal in J, Mj, and p but not 
diagonal in j and I. A similar transformation is discussed by 
Napolitano et al. l33ll in the context of ultra-cold collisions 
between atoms in the 1 Sq and 1 P 1 states. 

The next step is to express the body-fixed electronic 
states \SLjfl; a) = J2^. A (SLEA\jn}\SLi:A; a) in terms of 
\SLY,A; a), where S and A are projections of S and L along 
the internuclear axis. In other words, 

(S'L'j'n';a'\V iat \SLjn;a) = 

(S'L'Y,'A'\j'n'){SLT,A\jQ) 

E',A',£,A 

(S'L'E'A'; a'\V int \SLJ:A; a) . (A3) 

The Born-Oppenheimer potentials are diagonal in this 
\SLT,A;a) basis and the diagonal matrix elements are 
{SL^A;a\V int \SL^A;a) = VW( 2S+1 |A| CT ). 

Finally, we calculate matrix elements of the Hamiltonian 
Ha for two non-interacting atoms, where one atom is always 
in the 1 So state while the other atom can be in the state 1 So, 
1 Pi, or 3 P2,i,o- Because one atom is in the 1 So state, the 
molecular angular momenta S, L, and j in Eq. IAH are equiv- 
alent to those of the second atom. 
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TABLE I: The values of energies and lifetimes of Ca atomic states 
used in this paper. 



2s+ % 




"%■) [cm- 1 ] t a 


: 2S+ %o 


'Pi 




23 652.304 


4.59 ns 






15 315.943 




3 A 




15 210.063 


0.48 ms 


3 Po 




15 157.901 




'S 




0.000 





A realistic description of the atomic Hamiltonian should in- 
clude relativistic coupling between singlet 1 P 1 and triplet 3 Pi 
states. Therefore, the atomic Hamiltonian is not diagonal in 
the basis of Eq. JAli . Following Mies et a/. l39ll "dressed" 
electronic states for j = 1, which are a mixture of singlet and 
triplet, are introduced as follows 

\S = 0, Ljrrij; a) = 

cos(a)|5 = 0, Ljrrij; a) + sin(a) | S = l,Ljrrij;a) 

and 

| S = 1, Ljrrij; a) = 

— sm(a)\S = 0, Ljrrij; a) + cos(a)\S = 1, Ljrrij; a) 

where a is a small mixing angle. For j = and 2 we have 
\S Ljrrij; a) = \S Ljrrij; c). We assume that in this dressed 
basis the atomic Hamiltonian is diagonal with diagonal matrix 

elements (SLjlJMj;p\H A \SLjlJMj;p) = E A ( 2 ^ +1 Lj), 
where E A ( 2S+1 Lj) is the energy of the dressed state 2S+1 Lj 
relative to the 'So ground state. 

The mixing angle a is determined by the requirement 
that the dressed basis \S Ljrrij; a) reproduces the exper- 
imental transition probabilities between the excited 1,3 Pi 
states and the ground 1 So state. The angle can then 
be related to the ratio of the experimental dipole mo- 
ments of these transitions. In fact, we used tan(a) = 
^/{E a ^P 1 )/E a { 3 P 1 )] 3 [t a ^P 1 )/t a { 3 P 1 )}. This approach 
is only an appr oximation of the description of a real alkaline- 
earth atom ll39il . For the purposes of this paper, however, it 
is sufficient. The values for the e nerg ies and lifetimes have 
been obtained from Refs. 15311 and listed in Table |I] 

The experimental results on the 3 Pi lifetime are compared by 
Drozdowski et al. l53ll and some uncert ainty still exist. We 
have chosen the estimated value of Ref. IU9I1 . which lies be- 
tween the experimental values. 

The molecular Hamiltonian H conserves J, Mj, and p. 
In fact, the matrix elements of H in the basis of Eq. JA1I 
for given J and p are independent of Mj. It is con- 
venient to introduce the channels I7) = \SLjlJMj;p), 
I7) = \SLjlJMj;p), and note that ^Hff) = 
Y, 1 ', 1 {l'\l') l \l'\ H \l){l\l)- Close-coupling equations for 
the molecular wave function = \ j) F^{r)/r can be 



TABLE II: The allowed J e = 1, p R = -1 and M e = 0, ±1 
dressed channels of interacting ground- and excited-state alkaline- 
earth atoms. 



\%) 




L e 


je 


le 


1 


1 


1 


1 





2 


1 


1 


1 


2 


3 


1 


1 


2 


2 


4 





1 


1 





5 





1 


1 


2 



written as 

~Tu& ^M+E(7|ffA+Mnt+K.ot|7' > F 7'M = E F^r) , 
7' 

(A4) 

where E is the total molecular energy. These coupled 
Schrodinger equations are solved numerically to find scatter- 
ing and bound states. 

The collision between two ground state atoms can be solved 
separately from that of a ground plus an excited state atom. 
Moreover, for photoassociation of Ca we need scattering so- 
lutions of the Schrodinger equation for ground state collisions 
and bound states for ground-excited molecules. 

For the collision of two ground-state atoms 1 So+ 1 So there 
is only a single channel for given J g , M g and p g . The chan- 
nel is \-f g ) = \SgL g jgl g J g M g ;p g ) with S g = 0, L g = 0, 
j g = 0andZ g = J g , The two 40 Ca atoms are indistinguishable 
bosons and, therefore, the wavefunction must be symmetric 
under exchange of atoms and only channels with positive par- 
ity p = +1 exist. Consequently, only even l g = 0, 2, 4, . . . are 
allowed and the atoms interact on the gerade electronic state 
Solving the Schrodinger equation with E = e r > 
we obtain the scattering wavefunction |^+(£ r , J g , M g )) = 
\j g )F+ (r)/r. This wavefunction occurs in the stimulated 
width defined by Eq. (|5jl. For large r the wavefunction goes 

to Ff„(r) -»■ \/ 2/j 1 /(irh 2 k r ) sm(k r r + Trl g /2 + rn g ) exp(ir) lg ) 
IU4lHal . Our interest will be in the solution for J g = l g = 
or ,s-wave collisions. For collisions at ultralow temperatures 
other partial waves do not contribute significantly. 

For molecules formed by a ground- and excited-state atom 
there are multiple channels involved. The number of channels 
is determined by J e and p e . Photon selection rules limit the al- 
lowed total angular momentum of the excited bound states to 
Je = Jg, Jgil an d their parity to p e = —p g - For J g = or s- 
wave collisions PA can only make J e = 1 andp e = — 1 bound 
states. Tablelllllists the five channels for J e = 1 andp e = —1. 
The first two channels correspond to 1 So+ 3 Pi states. There 
are no channels with j e = and, since the atomic energies 
satisfy E A ( 3 P ) < E A ( 3 Pi) < E A ( 3 P 2 ), predissociation of 
bound states below the 1 So + 3 Pi limit does not occur. In fact, 
this is true for all odd J e and p e = —1. 

The numerical solutions of the Schrodinger equation for 
the scattering collision in the grou nd state are calculated us- 
ing the Numerov method 1541 15511 implemented in the close 
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coupling code developed by Mies, Julienne and Sando |5( 
The coupled-channel bound states calculations for the excited 
state are carried out using the discrete variable representation 
(DVR) HH following Tiesinga et al. 01- 



APPENDIX B: INTERACTION WITH LIGHT 

In our treatment of the photoassociation process we need 
the stimulated width, defined by Eq. ®, between the scatter- 
ing ground state and the excited bound states. The stimulated 
width can be expressed as 



2tt 



E 

7e 



dr (%\V las \ lg )F~*(r)F+(r) 



(Bl) 

where the wavefunctions are expressed in terms of channel 
functions F^ e (r) and F+(r) and (7 e |V!as|7g) are matrrx el- 
ements between channels in the ground and excited state. If 
we only consider dipole transitions and neglect retardation ef- 
fects, this matrix element is independent of the interatomic 
separation r and given by 



(7e|Vla S |7s) 



2nl 1 

C V47T£0 



(%\d-e q \ lg ) (B2) 



where the laser has intensity / and polarization e q , eo the per- 
mittivity of vacuum, and (%\d- e q \-f g ) = ^-yjjehe) (le\d ■ 
e q \-f g ) is the molecular dipole matrix element between spin 
channels. 

Napolitano et al 13311 discuss the connection between 
molecular and atomic dipole matrix elements. A similar ap- 
proach is used here to find matrix elements. The molecular 
dipole matrix elements between ground and excited state are 

(S e L e j e l e J e M e \p e \d q < \SgL g j g l g J g M g ;p g ) = 



1 (2/ e + l)(2l g + 1) 



(J g KM g q\J e M e 



(2J e + l)2 

]T (j g Kn g q'\j e n e )(j e i e n e o\j e Q e )(j g i g Q g o\j g n g ) 

(S e L e j e n e ;a e \d£, \S g L g j g n g :a g ) (B3) 

where the spherical tensor operator d3 = d ■ e q with K = 1 
and q = —1,0, or +1. In Eq. JB3I > we have expressed 
the dipole operator in a body-fixed coordinate system. For 
both singlet S — ground and excited states we realize that 
\SLjCt; a) = \SLT,A; a) and thus j = L and fl = A. It then 
follows 

lK 



(S e L e j e Q e ; o e \d q \S g L g j g Q g ; a g 



'I 

(S e L e Z e A e ;a e \df\S g L g Z g A g ;cr g ) , 



for electronic states with zero total electron spin Sq - — ■ Sg — 0. 
In the \SLY,A; a) basis the dipole operator can be evaluated 
in terms of the atomic line width, TaQ~Pi), of the 1 P\ atomic 



state decaying to the ground 1 Sq state. There are four distinct 
matrix elements (010A; a\d\\0000; g) = d^A^) given by 



dC\M<r) = 




4 7 re r M ( 1 |A| (T ), (B5) 



where rA/( 1 |A| (T ) is the molecular linewidth. The d ^lAL, 
are r-independent quantities as we neglect retardation 1191148 
and, therefore, the molecular linewidth can be well ap- 
proximated by r^/( 1 |A| ll ) = 2Ta( 1 Pi) for ungerade states 
and rM( 1 |A| s ) = for gerade states. In our phase conven- 
tion the <i( 1 |A| cr ) are positive. The molecular dipole matrix 
elements between the triplet S e = 1 and the ground S g = 
states are zero. 



APPENDIX C: NEAR THRESHOLD VIBRATIONAL 
QUANTUM NUMBER 



Following LeRoy and Bernstein 1 60] the JWKB quan- 
tum condition for the eigenvalues of energies E for a po- 
tential V(r) is that Q(E) = vir where v is an inte- 
ger and the phase <&(E) = Qjwkb(E) - n/2. Here 
Qjwkb(E) = dr k(E, r) is the semiclassical phase in- 
tegral of k{E, r) = \J (E — V(r))2fi/h calculated between 
the classical turning points r\ and r 2 . This criterion is valid 
for deepl y bo und states but breaks down near the dissocia- 
tion limit I61B . It can be shown that for energies E close to 
the threshold this condition still holds l62Tl if the phase <f>(E) 
is modified in the following way: = $>jwkb(E) — 

n/2 — n/(2q — 4) for long-range potentials of the form 
V(r) = C q /r q . A recent discussion on improvements to the 
LeRoy-Bernstein approach can be found in Ref. l63ll . 

The phase $>(E) allows us to define a generalized vibra- 
tional quantum number as ($(£') — Vq^/tt, a continuous 
function of energy E, where vo is an arbitrary constant, which 
does not need to be an integer. For vq = the generalized vi- 
brational quantum number is 0, 1, ... ,v max — 1, or iJ max when 
E is equal to the energy of a bound state. Here, a zero value 
corresponds to the most deeply bound state and v ma _ K is the vi- 
brational quantum number of the last bound state. At the dis- 
sociation limit this generalized vibrational quantum number 
equals vp = $(0)/ir and, in general, can have a non-integer 
value II60B. 

In the analysis of bound states near the threshold it is con- 
venient to set vo = vu or Vq = v max + 1. In the first case 
the generalized vibrational quantum number is equal to zero 
at threshold. Negative integers -1, -2, -3, ... correspond to en- 
ergies that mark the edges of "bins", in which there is exactly 
one bound state. Only for integer vjj, however, "bin edges" 
coincide with bound state energies of the potential. In fact, 
"bin edges" were introduced in this way in Fig.|4] A second 
useful choice of Vq is vo = v max + 1. It allows us to de- 
fine a near threshold vibrational quantum number, which has 
values -1, -2, -3, ... for bound states counting from the top. 
The near threshold vibrational quantum number for E = is 
vd — v ni ax — 1> and lies between -1 and 0. 
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-10 o 

A e /h [GHz] 

FIG. 10: The near threshold vibrational quantum number u(A e )/n 
of the potential V( 3 I1 U ) + [2 + J e (J e + l)]h 2 /(2fj,r 2 ) with J e = 1 
as a function of the binding energy A e . The dashed lines indicate the 
binding energies at which ^(A e )/7r is integer. The value of v(A e )/n 
at A e = is also indicated. 



In practice, for energies far from threshold but not far 
enough to justify the use of the JWKB approximation the en- 
ergy dependence of the effective vibrational quantum num- 
ber can be calculated using expressions given by Mies Eh . 
Following Refs. Ir34l l65ll we define a near threshold vibra- 
tional quantum number i/(E)/n as a continuous function of 
energy E, where v(E) is the phase difference between two 
solutions of the Schrodinger equation /(r) and <p(r) at the 
equilibrium separation, r e , of the potential V(r). The func- 
tion f(r) is obtained by solving the Schrodinger equation as- 
suming that /(0) = 0. The function <j>(r) is obtained by solv- 
ing the Schrodinger equation assu ming that <f>(oc) = 0. Then 
tan(f(.E)) can be calculated from 16411 



tan(v(E)) 



K(E,r e )[f(r e )<j>\r e )-<l>(r e )f'(r e )\ 
^iE,r e )f(r e )<t>ire) + <t>'(r e )f'(r e ) 



(CI) 



where the primes denote the first derivative with respect to r. 
The near threshold vibrational quantum number, defined this 
way, has an integer value for energies corresponding to the 
bound states of the potential V(r). 

A connection between the near threshold vibrational quan- 
tum number and the binding energies A e of the multichannel 
excited bound states, used in Fig. [8] can be made using quan- 
tum defect theory 16411 . The theory states that the shape of the 
energy dependence of v{E) near threshold is nearly indepen- 
dent of the short-range form of the potential and that we can 
replace the energy E in Ea. lCll bv the binding energy A e . 

The near threshold vibrational quantum number ^(A e )/7r 



for a potential F( 3 n u ) + [2 + J e (J e + l)]?i 2 /(2^r 2 ) as a 
function of A e is shown in Fig.^3 F° r this potential the short 
range shape is adjusted to obtain single-channel bound states 
(or negative integer near threshold vibrational quantum num- 
bers) for binding energies equal to -1.39 GHz, -8.71 GHz, - 
25.56 GHz and -54.01 GHz. These values agree well with 
the multi-channel "bin" edges -1.45 GHz, -8.75 GHz, -25.56 



1000r 



1 2 3 

Icxrjr ( c )/r J 
. 10 eg X o.rat' 




FIG. 1 1 : The stimulated linewidth of the 0„ band as a function of the 
near threshold vibrational quantum number vi A e )/n and the scatter- 
ing length a sca t of the ground state potential. The laser intensity is 1 
W/cm 2 and the kinetic energy is e r /fcs = 1 AtK. 



GHz and -53.94 GHz shown in Fig.|4] In fact, for binding en- 
ergies below -0.026 GHz the multi-channel 0+ bound states 
are well approximated by this single-channel Hund's case(c) 
potential. The connection between integer i/(A e )/7r and the 
multi-channel bin edges breaks down for smaller binding en- 
ergies as Coriolis mixing changes the multi-channel coupling 
scheme from Hund's case (c) to (e). Since a single-channel 
can not fully simulate the multi-channel bound states, the near 
threshold vibrational quantum number ;y(A e )/7r = —0.2135 
is not integer for A e = 0. 

Figure ITTI shows the results of Fig. [8] in terms of the near 
threshold vibrational quantum number v{ A e )/7r instead of the 
binding energy A e . It is clearly seen from Fig. ^2 that f° r 
j/(A e )/7r between v and v + 1, where v = —4, —3, or — 2, 
and thus for each corresponding bin a pair of dark structures 
occurs. In fact, the "first kind" of structure always appears at 
nearly the same location within a bin. The "second kind" of 
structure moves much more within a bin. Both observations 
reflect the fact that the Cq coefficients of the ground and exited 
state potentials are similar. 
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